Non-resonant driving of H atom with broken time-reversal symmetry 
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The dynamics of atomic hydrogen placed in a static electric field and illuminated by elliptically 
polarized microwaves is studied in the range of small field amplitudes where perturbation calculations 
are applicable. For a general configuration of the fields any generalized time-reversal symmetry is 
broken and, as the classical dynamics is chaotic, the level statistics obeys the random matrices 
prediction of Gaussian unitary ensemble. 
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The Bohigas, Giannoni and Schmit conjecture Q that 
a general quantum system with underlying chaotic clas- 
sical dynamics has statistical properties of energy levels 
described by random matrices theory has been confirmed 
by theoretical studies of numerous chaotic systems, see 
e.g. [||]|]. Experimental verifications on the other hand 
are much less abundant and, as far as we know, 

concern only systems possessing anti-unitary (general- 
ized time-reversal) symmetry where the level statistics 
can be modeled by random matrices of the Gaussian or- 
thogonal ensemble (GOE) ||. Realization of a quantum 
system with quadratic level repulsion as is typical for 
matrices from the Gaussian unitary ensemble (GUE) ||] 
requires breaking of any anti-unitary symmetry invari- 
ance. E.g. for atomic or molecular systems it means 
one would have to apply magnetic field inhomogeneous 
across the molecule being well experimentally controlled 
on such a small scale Q. Such a requirement is rather 
unattainable even in the Rydberg regime of excitation. 
While a real quantum system with GUE statistics has not 
been realized experimentally, there are experiments with 
microwave cavities (so-called wave chaos experiments) 
where anti-unitary symmetry can be broken by apply- 
ing some ferrite devices |7],|8| . 

Recently, however, it has been shown |j|,[l0| that to 
break anti-unitary symmetry invariance, in atomic sys- 
tems, it is not necessary to employ magnetic field. Indeed 
a combination of elliptically polarized microwaves with a 
static electric field applied to, e.g., hydrogen (H) atoms 
can do this work as well. In the previous studies ||,[l0) we 
have restricted ourselves to a case when the microwave 
field is resonant with the classical motion of the electron, 
i.e. the ratio of microwave frequency, to, to Kepler fre- 
quency, lok — V n o (where no is the principal quantum 
number of H atom), is an integer number. It has been 
found that such a system can reflect level statistics very 
close to the GUE prediction. 

In this paper we study a more general situation of non- 
resonant driving of the atoms and show that appropriate 
choice of the system parameters allows us to reach sta- 
tistical properties expected for matrices of the GUE. We 
consider weak fields limit where quantum results may be 
obtained by the lowest orders of the perturbation theory. 



Classical investigation of the system behavior is carried 
out also in terms of the perturbation calculations. It 
allows us to find out for which system parameters the 
secular motion of the electronic ellipse reveals chaotic 
dynamics. 

The Hamiltonian of a realistic three-dimensional H 
atom placed in a static electric field and driven by an 
elliptically polarized microwave field in atomic units, ne- 
glecting relativistic effects, assuming infinite mass of the 
nucleus, and employing dipole approximation reads: 
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F(xcosiot + ay sincot) + E • r, (1) 



where F, a and u) are, respectively, the amplitude, degree 
of elliptical polarization and frequency of the microwave 
field while E stands for the static electric field vector. 
As the external perturbation is time periodic, one may 
apply Floquet formalism and look for quasiencrgy levels 
of the system. 

We assume very small field amplitudes which allows in- 
cluding only the lowest non- vanishing terms in the quan- 
tum effective Hamiltonian ijllf describing the dynamics 
of a given uq hydrogenic manifold. The static electric 
field contributes already in the first order of E as the 
states with fixed no are coupled among each other by 
the E • r operator. On the other hand there is no direct 
coupling between the states due to the microwave pertur- 
bation. Indeed the states can be coupled only indirectly 
by the process of absorption and emission of microwave 
photons. So that the first non-vanishing term is second 
order in F . The final matrix (with n\ dimension) of the 
quantum effective Hamiltonian, i.e. the Hamiltonian in 
the first order in E and second order in F, is then diago- 
nalized by standard routines. The details of the quantum 
perturbation calculations can be found in [juj, here we 
would only like to stress that the whole physically real- 
istic problem is reduced to the analysis of the finite n§ 
dimensional Hilbert space where 1/no plays the role of 
the effective Planck constant. 

We are now turning to classical analysis of the system. 
The range of our interests is the high frequency regime, 
i.e. uj > lok, which means that for weak external fields 
we have two fast degrees of freedom in the system, i.e. 
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the position of the electron on an elliptical orbit and the 
phase of the microwave field, and two slow degrees of 
freedom corresponding to the orientation of the elliptical 
trajectory. One can get rid of the fast degrees of freedom 
by means of the classical perturbation theory [M which 
results in the classical effective Hamiltonian describing 
slow precession of an electronic ellipse. The perturbation 
calculation can be easily carried out employing the Lie 
method |12| , actually it closely follows the similar proce- 
dure applied in [|l3| to the H atom perturbed by linearly 
polarized microwave field. The first stage is to express 
the Hamiltonian (Q) in terms of the action-angle variables 
of the unperturbed hydrogen atom. The new pairs of the 
canonically conjugate variables are (J, 0), i.e. principal 
action (analog of the principal quantum number, no) and 
position of the electron on an ellipse respectively; (L, $), 
i.e. angular momentum of the electron and conjugate 
angle; and finally (Af, $), i.e. angular momentum pro- 
jection on z axis and angle of rotation around this axis 
J[4| |. Averaging the resulting Hamiltonian over 9 and 
t immediately gives the first order contribution to the 
classical effective Hamiltonian 



i/ (1) (L,*,M,$) 
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The second order contribution of microwave field requires 
calculating the generating function, w p2| which is the 
solution of the following equation 



dw 
~dt 



dw 
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where H m i cro is the microwave part of the Hamiltonian 
(|l|) expressed in the action-angle variables (explicit ex- 
pression for Hmicro as a Fourier series can be found in 
PI). The solution for w is given as an infinite series 
with terms containing l/(w ± mux), where m is an in- 
teger number. For resonant driving, i.e. uj/ujk ~ m, 
one faces the small denominators problem p2[ but here 
we are not affected by this problem as we are interested 
in a non-resonant perturbation. Having calculated w it 
is straightforward task to obtain the second order con- 
tribution to the effective Hamiltonian by averaging the 
Poisson bracket of w and H m i cro over and t, i.e. 



H^(L,^>,M,<i>) = ^({w,H mlcro }) 



e,i> 



(4) 



(we omit the lengthy explicit formula for here). The 
final classical effective Hamiltonian reads 



H e H (L, M, $) = ~+H {1) +H {2) . (5) 

This is the classical counterpart of the quantum effec- 
tive Hamiltonian, namely first order in the static electric 
field and second order in microwave field. The classical 
Hamiltonian (|s|) possesses scaling symmetry, i.e. one can 
get rid of one of the parameters of the system. Intro- 
ducing F a = nfjF, E = UqE, uj = tiqU, L = L/n a , 
Mq = M/na and = n 2 ) H e ff the dynamics becomes 

independent of the particular choice of the no hydrogenic 
manifold. 

For a general fields configuration the secular motion 
in the (L, ^, M, $) phase space is not integrable and to 
investigate classical dynamics we have to perform numer- 
ical integration of the equations of motion generated by 
the Hamiltonian (^|). For the linear microwave polariza- 
tion without additional static electric field considered in 
|l3"f the secular motion was one dimensional and employ- 
ing the WKB quantization rule the authors were able 
to get semiclassical predictions for quasienergy levels in 
a very good agreement with exact quantum numerical 
data. 

For elliptically polarized microwaves and general ori- 
entation of the static field vector the anti-unitary sym- 
metry invariance is broken. Only when either E x — 
or E y — the system is invariant with respect to the 
time-reversal combined with x — > — x or y — * —y trans- 
formations respectively, see ([!]). As an example of a 
general elliptical polarization case we have further on 
analyzed the degree of the polarization a — 0.4. The 
microwave frequency has been chosen as luq = 1.304 
and the amplitude as Fq = 0.02 which is well in the 
range where, for the linear microwave polarization, the 
WKB calculations give good agreement with the exact 
numerical results |l3| . Then by investigating Poincare 
surface of section we have found that the amplitude 
Eq = 0.00028 and the orientation of the static field vec- 
tor (f> ks 0.37T, 9 w 7r/4 (where 9, (j) are usual spheri- 
cal angles) correspond to chaotic dynamics in the energy 
interval H^ ff e [-0.50032,-0.50008]. Putting = 
(E y = 0) one recovers anti-unitary symmetry of the sys- 
tem, then the classical dynamics is found to be predom- 



inately chaotic for H. 
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-0.50016,-0.49992]. Fig. | 
shows examples of the phase space structures for both 
= and (f) — 0.37T. 

Having done classical analysis we can switch to quan- 
tum calculation results. In order to get reasonable statis- 
tics for quasienergy levels we have diagonalized the ma- 
trix of the quantum effective Hamiltonian for different 
no manifolds in the range no = 50 -j- 59. From each di- 
agonalization we have separated and unfolded spectrum 
in the energy intervals corresponding to chaotic classi- 
cal dynamics. This procedure has been applied to the 
anti-unitary invariant case (<f> = 0) and to the case with 
broken anti-unitary symmetry (</> = 0.37r). Fig. |^ shows 
histograms of the nearest neighbor spacing (NNS) dis- 
tributions of the quasienergy levels (there are about 10 4 
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spacings in each of the data sets) together with the plots 
of the Wigner surmises § for the GOE and GUE. The 
figure also shows the spectral rigidities, i.e. A3 statistics 
H . The qualitative agreement of the data with the ran- 
dom matrices theory is apparent, especially for the case 
with broken anti-unitary symmetry. 

To focus on quantitative measure we have fitted the- 
oretical NNS distributions to the data (to avoid the de- 
pendence of the results on bin size the distributions have 
been fitted to the cumulative histograms ]l5|] ). For the 
anti-unitary invariant case the best fitting Berry-Robnik 
distribution (i.e. the distribution for independent super- 
position of Poisson and GOE spectra jl6| ) results in the 
parameter value (relative measure of the chaotic part of 
the phase space) q = 0.98 with x 2 /N — 0.3, i.e. chi- 
squares divided by the number of spacings N. Employing 
the fzrailev distribution Jf7];[l8| we get the value of the 
levels repulsion parameter f3 = 0.96 with % 2 /N = 0.8. 
For the spectral rigidity the best fitting A3 statistics 
corresponding to an independent superposition of Pois- 
son and GOE spectra |3| results in q = 0.99 which is 
in agreement with the value obtained from the Berry- 
Robnik distribution fit. 

In the case with broken anti-unitary symmetry one gets 
the following parameters values: Berry-Robnik statistics 
(now it corresponds to independent superposition of the 
Poisson and GUE spectra ||) q = 1 with X 2 /N = 0.2, 
Fzrailev distribution (3 ~ 2.05 with x 2 /N = 0.3 and from 
the spectral rigidity fit q = 1. 

The presented results show undoubtedly that the sys- 
tem under consideration can reflect the Bohigas, Gian- 
noni and Schmit conjecture both with or without anti- 
unitary symmetry. The question is if such a behavior can 
be observed experimentally? With this respect let us 
consider the case with broken generalized time-reversal 
symmetry As there is no discrete symmetry in the sys- 
tem one has no problems with separation of overlapping 
spectra. On the other hand it results in a big density 
of states and requires high experimental resolution. We 
have presented the data for hydrogenic manifolds in the 
range no = 50 -7- 59 in order to have good statistics but 
our results for no = 40 -f- 49 reveal the very same be- 
havior. For example for no — 40 and F — 40 V/cm, 
w = 2tt • 134 GHz, E = 0.56 V/cm the average level 
spacing is of order 10~ 4 crn^ 1 and a measurement in the 
energy range between —68.628 cm -1 and —68.595 cm -1 
provides about 500 levels. It sounds feasible experimen- 
tally. 

The previous studies [p|Jl0[ have been devoted to res- 
onant microwave driving of H atoms placed in a static 
electric field. The present work deals with a general non- 
resonant driving and shows that intra-manifold chaotic 
dynamics ]2C|-|2^J^,p^| , i.e. the situation when states cor- 
responding to a fixed principal quantum number no are 
mixed significantly only among each other and under- 
lying classical dynamics is irregular, is not restricted to 
a particular resonant case but exists widely in the fre- 
quency domain. The presented behavior is field ampli- 



tudes independent, i.e. if F — > and E — > but E/F 2 — 
const the structure of the phase space corresponding to 
secular motion remains unchanged as are the statistical 
properties of quasienergy levels. This is a signature of 
inapplicability of the Kolmogorov-Arnold-Moser theorem 
|l2| to a highly degenerate Coulomb problem. For high 
field amplitudes when hydrogenic manifolds can not be 
considered isolated and the inter-manifold mixing comes 
to the scene our perturbation approach, obviously, be- 
comes irrelevant. However this regime, for a realistic 
three-dimensional problem, still constitutes a challenge 
for the theory. 

The author is grateful to Dominique Delande and 
Jakub Zakrzewski for discussions. The results were at- 
tained with the assistance of the Alexander von Hum- 
boldt Foundation. The partial support of KBN under 
project 2P302B00915 is acknowledged. 



[1] O. Bohigas, M. J. Giannoni, and C. Schmit, Phys. Rev. 
Lett. 52 (1984), 1. 

[2] F. Haake, Quantum Signatures of Chaos, Vol. 54 of 
Springer Series in Synergetics Springer, Berlin, 1991. 

[3] O. Bohigas in Chaos and Quantum Physics, 1899 Les 
Houches Lecture Session LII, edited by M.-J. Giannoni, 
A. Voros, and J. Zinn- Justin (North-Holland, Amster- 
dam, 1991), p. 87. 

[4] T. Zimmermann et ai, Phys. Rev. Lett. 61 (1988), 3. 

[5] H. Held, J. Schlichter, G. Raithel, and H. Walther, Eu- 
rophys. Lett., 43, 392 (1998). 

[6] M. L. Mehta, Random Matrices Academic Press, San 
Diego, 1991. 

[7] P. So, S. M. Anlage, E. Ott, and R. N. Oerter, Phys. Rev. 

Lett. 74 (1995), 2662. 
[8] U. Stoffregen, J. Stein, H. J. Stockmann, M. Kus, and F. 

Haake, Phys. Rev. Lett. 74 (1995), 2666. 
[9] K. Sacha, J. Zakrzewski, and D. Delande, Phys. Rev. 
Lett. 83 (1999), 2922 and 84, 200 (2000). 

[10] K. Sacha, J. Zakrzewski, and D. Delande, Ann. Phys. 
(N.Y.), in press. 

[11] C. C. Cohen- Tannoudji, J. Dupont-Roc, and G. Gryn- 
berg, Atom-Photon Interactions: Basic Processes and 
Applications John Wiley and Sons, New York, 1992. 

[12] A. J. Lichtenberg and M. A. Lieberman, Regular and 
Stochastic Motion, Vol. 38 of Applied Mathematical Sci- 
ences Springer, Berlin, 1983. 

[13] A. Buchleitner and D. Delande, Phys. Rev. A 55, R1585 
(1997). 

[14] K. Sacha and J. Zakrzewski, Phys. Rev. A 58 (1998), 
488. 

[15] T. Prozen and M. Robnik, J. Phys. A 26 (1993), 2371. 
[16] M. V. Berry and M. Robnik, J. Phys. A 17, 2413 (1984). 
[17] F. M. Izrailev, Phys. Rep. 196 (1990), 299. 
[18] G. Casati, F. M. Izrailev, and L. Molinari, J. Phys. A 24 
(1991), 4755. 



3 



[19] M. Robnik, J. Dobnikar, and T. Prosen, J. Phys. A 32, 
1427 (1999). 

[20] J. von Milczewski, G.H.F. Diercksen and T. Uzer, Phys. 

Rev. Lett. 73 (1994), 2428. 
[21] J. von Milczewski and T. Uzer, Phys. Rev. E 55 (1997), 

6540. 

[22] J. Main, M. Schwacke, and G. Wunner, Phys. Rev. A 57 
(1998), 1149. 




.... . ■ 



I 



FIG. 1. Poincare surface of section (at $ = 0) of the clas- 
sical secular motion, Eq. (^), for the hydrogen atom placed in 
a static electric field, with the amplitude Eo = 0.00028, and 
driven by an elliptically polarized (a = 0.4) microwave field 
with frequency u)q = 1.304 and amplitude Fo = 0.02. The 
coordinates used for the plot are the scaled angular momen- 
tum Lo = L/no and its canonically conjugate angle ^. Panel 
(a) corresponds to the anti-unitary invariant case, i.e. the 
orientation of static field vector is <j> = 0, if) = n/4; with the 
energy Hq = —0.50008. Panel (b) is related to the broken 
anti-unitary invariance, i.e. <f) = 0.3-7T, ifr = 7r/4; with the en- 



ergy H 



-0.5002. Note that, for the parameters chosen, 




FIG. 2. Nearest neighbor spacing distribution and spec- 
tral rigidity, A3 statistics, for hydrogen atom placed in a 
static electric field and illuminated by microwave field, for 
the case with [panels (a)-(b)] and without [panels (c)-(d)] 
anti-unitary symmetry invariance, compared with predictions 
of random matrices ensembles. The data have been col- 
lected for no = 50 4- 59 with the same fields parameters as 
in Fig. ^. In panels (a) and (c) solid lines indicate the best 
fitting Izrailev distributions, while dashed and dash-dotted 
lines correspond to GOE and GUE distributions respectively 
(in panel (c) the dash-dotted line is invisible behind the solid 
one). Panels (b) and (d): solid and dotted (hardly visible be- 
hind the solid lines) lines correspond to numerical data and 
their best fits, while dashed lines indicate Poisson, GOE and 
GUE predictions as indicated in the figure. 



not the whole (Lo,ip) space is accessible. 
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